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Modelling fully convective stars in eclipsing binaries: 
KOI-126 and CM Draconis 
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ABSTRACT 

We present models of the components of the systems KOI-126 and CM Draconis, 
the two eclipsing binary systems known to date to contain stars with masses low 
enough to have fully convective interiors. We are able to model satisfactorily the 
system KOI-126, finding consistent solutions for the radii and surface temperatures 
of all three components, using a solar-like value of the mixing-length parameter a in 
the convection zone, and PHOENIX NextGen ID model atmospheres for the surface 
boundary conditions. Depending on the chemical composition, we estimate the age 
of the system to be in the range 3-5 Gyr. For CM Draconis, on the other hand, wc 
cannot reconcile our models with the observed radii and T c g using the current metal- 
poor composition estimate based on kinematics. Higher metallicities lessen but do 
not remove the discrepancy. We then explore the effect of varying the mixing length 
parameter a. As previously noted in the literature, a reduced a can be used as a simple 
measure of the lower convective efficiency due to rotation and induced magnetic fields. 
Our models show a sensitivity to a (for a < 1.0) sufficient to partially account for 
the radius discrepancies. It is, however, impossible to reconcile the models with the 
observations on the basis of the effect of the reduced a alone. We therefore suggest that 
the combined effects of high metallicity and a reduction could explain the observations 
of CM Draconis. For example, increasing the metallicity of the system towards super- 
solar values (i.e. Z = 2 Zq) yields an agreement within 2 a with a = 1.0. 
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1 INTRODUCTION 

Low-mass stars (i.e., M < 1 Mq) represent the largest frac- 
tion of the stellar population in the Galaxy. The availability 
of a reliable theoretical mass-radius relationship is of con- 
siderable interest both per se and in view of the applica- 
tion to the study of exoplanet-hosting stars. Moreover, be- 
low about 0.30 Mq, the depth of the subsurface convection 
zone extends enough to reach the centre and stars become 
fully convective (hereafter FC). This new regime has pe- 
culiar characteristics and presents additional challenges. In 
constructing theoretical models of FC stars, special care is 
required in handling all the main input physics ingredients, 
such as the equation of state, the opacities, and the atmo- 
spheric boundary conditions, due to the lower temperature 
and h ighe r density regime (see the reviews bv lAllard et al.l 
Il997l and IChabrier fc Baraffel |200Ch . Although convection 
becomes very efficient at quite shallow depths (i.e. V ~ V a d 
at T > T c ff, see e.g. section 7.3.3 of [Hansen et al.ll2004 ). 
it is st ill necessary to rely on the mixing length theory 
(MLT, iBohm-Vitensj 1 19581 ) for the convective transport 



in the outermost layers and in the atmosphere, which can 
have a significant impact o n the global structure parameters 
|Chabrier fc Baraffdll997h . 

Double-lined eclipsing binary systems (hereafter EBs) 
provide very accurate determination of the fundamental pa- 
rameters of stars, allowing to test theoretical model predic- 
tions. There is mounting evidence that low-mass stars in 
EBs have radii about 10 per cent larger and effective tem- 
per atures about 5 per cent cooler than theoretical mo dels 
(see iRibafeod ; iRibas et al]|20Qg| ; iMorales et all 12009] and 
references therein). Interestingly, the discrepancies are in 
the right proportions to pro duce roughly the sam e luminos- 
ity. To explain these results. IChabrier et al.l (|2007l ) proposed 
that close EBs may not be representative of the whole low- 
mass stars population because of the high rotation regimes 
produced by spin-orbit synchronization via tidal interaction. 
This explanat ion is corroborated by the significant correla- 
tion found by lLopez- Morales! (|2007l ) between the difference 
in the radii and the level of magnetic activity. Indeed, both 
rotation and magnetic fields can reduce the efficiency of con- 
vection, which is measured by the parameter a in the stan- 
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dard MLT framework (e.g. lTavlerlll973h . A faster rotation, 
moreover, can enhance the magnetic activity through the 
dynamo effect, leading to a substantial increase in the star- 
spots coverage of the surface. Both effects were phen omeno- 
logically taken into account bv lChabrier et al.1 (|2007T l. using 
a lower a with respect to the solar-calibrated value and re- 
ducing the radiative flux at the surface proportionally to the 
star-spots filling factor /3, respectively. Their results show 
that the parameters a and /3 have a degenerate effect on 
the stellar radius. They also claim that FC stars are almost 
insensitive to a variation of a and very low values (e.g., 
a < 0.1) are required to reconcil e the models with th e ob- 
servations. Finally, according to iMorales et ail (|2010h . the 
discrepancy may be partially accounted for by a systematic 
error in the light curve-derived radii, which would occur if 
a value of /3 of about 35 per cent and a conce ntration of the 
spots near the stellar poles are postulated. As lBurrows et ail 
(|201ll ) have pointed out, however, theoretical models calcu- 
lated with a quite restricted set of parameters (e.g., solar 
metallicity), have been compared with observations of stars 
of poorly determined or even unknown metallicity, assum- 
ing a general validity that is not rigorously justified. Even 
a moderate increase in metallicity can result in theoretical 
radii increases which are significant in comparison with the 
discrepancies discussed so far. A combination of all these ef- 
fects is thus the most likely solution to the low-mass stellar 
radii problem. 

In this paper, we attempt to model the components of 
the two EBs known to date to host FC stars. Our main focus 
is to determine the sensitivity of the best-fitting models to 
various possible choices of the parameters, with as much gen- 
erality as possible. The paper is organized as follows: in Sec- 
tion[5]we discuss the available data and compare the physical 
properties of the two systems; in Section [3] we present the 
code and input physics used in our calculations and in Sec- 
tion 13.11 we show the impact of the atmospheric boundary 
conditions on our models. In Sections 14.11 and 14.21 we dis- 
cuss our best-fitting solutions for the systems KOI- 126 and 
CM Draconis (hereafter CM Dra), respectively. We compare 
our re sults for KOI-126 with those obtained bv lFeiden et al.l 
|201ll) in Section[473] Finally, in Section[5j we review results 
and conclusions of this work. 



2 THE SYSTEMS KOI-126 AND CM DRA 

The systems KOI-126 (KIC number: 5897826) and CM Dra 
host the only EBs known to date in the FC mass range. 
The former was re cently discovered by the Kepler mission 
(|Carter et al.|[201ll ') ; it is a hierarchical triple system, with a 
low-mass binary subsystem (orbital period: 1.767 d) orbiting 
the primary sta r with a period of 33. 92 d. The second system 
(e.g. lLacvlll977l . IMorales et al.ll2~009l . and references therein) 
is composed of two FC stars (orbital period: 1.268 d) which 
have common proper motion with a white dwarf. 

For the stars in these systems, masses and radii are 
determined with very high precision from the modelling of 
the light curve eclipses. Other parameters however, such as 
metallicity and effective temperatures, are much less well 
known. Table [T] summarises the known physical parameters 
of the two systems. 

Spectroscopic observations of KOI-126 only show fea- 



Table 1. Physical parameters of the two systems. The KOI-126 
data marked with a t are estimates based on the quoted flux 
ratio of each component to the primary: /b//a = 3.26 ± 0.24 ■ 
10~ 4 , Self A = 2.24 ±0.48 ■ 10~ 4 (sec text). The orbital periods 
are: P A BC = 33.92 d (KOI 126 B-C binary around KOI-126 A); 
P BC = 1.767 d (KOI-126 B-C binary); P 12 = 1.268 d (CM Dra). 
The metallicity range of CM Dra reported here is that given in 
the literature. 





KOI-126 




CM Dra 




M 


ass [Mq] 




M A 


1.347 ± 0.032 






M b 


0.2413 ± 0.0030 


Mi 


0.2310 ± 0.0009 


M c 


0.2127 ± 0.0026 


M 2 


0.2141 ± 0.0010 




Radius [Rq] 




Ra 


2.0254 ± 0.0098 






Rb 


0.2543 ± 0.0014 


Ri 


0.2534 ± 0.0019 


Rc 


0.2318 ± 0.0013 


R2 


0.2396 ± 0.0015 




Log. surface gravity, log 10 


(g/cm 2 s _1 ) 


log<U 


3.9547 ± 0.0069 






l°g 9B 


5.0101 ± 0.0029 


log gi 


4.994 ± 0.007 


log gc 


5.0358 ± 0.0027 


log S2 


5.009 ± 0.006 




Effective temperature [K] 


T e B,A 


5875 ± 100 






T e g,B 


(3300)+ 




3130 ± 70 


Teff,C 


(3200)+ 


T e ff,2 


3120 ± 70 




Luminosity 






0.6417 ± 0.029 






(-2.14) + 
(-2.26) + 




-2.258 ± 0.038 
-2.313 ± 0.056 




Metallicity 




[Fe/H] 


+0.15 ± 0.08 


[M/H] 


-1.0 < [M/H] < -0.6 



tures associated with the primary star (|Carter et alJlioill . 
see also the online supporting material). However, the light 
curve analysis provides a determination of the fluxes of 
the two low-mass stars relative to the primary, which can 
be used to derive a rough estimate of their effective tem- 
peratures. We proceeded as follows: using synthetic spec- 
tra from the NextGe n (NG) grio3 (construc ted with the 
PHOENIX ID code, iHauschildt et all I1999T ), with solar 
metallicity and appropriate logg and logT c ff, we calculated 
the theoretical integrated flux for each star of the sys- 
tem, taking into account the response function of Kepler 
(|Van Cleve fc CaldwelT1l2009l ). We estimate that the effec- 
tive temperatures fall in the range 3200-3300 K; since the 
radii are known, the luminosities can be derived as well (see 
Tabled]). 

Although considerable effort has been devoted to the 
study of CM Dra, both spectroscopic ally and photometri- 
cally (|Lacvlll977l . iMetcalfe et al.lll996l ). the direct determi- 
nation of its metallicity is still controversial l|Viti et al ] |l997i . 
120021 : IMorales et alll2009l ): the system is considered a Pop- 
ulation II member based on its measured proper motion. 

The masses of the FC components of both systems fall 
in the range 0.20 < M/M@ < 0.25 and have remarkably 
similar ratios. Interestingly, KOI-126 C and CM Dra 2 have 
the same mass within the error, but the radius of the latter is 
significantly larger (of about 3 per cent). Similarly, KOI-126 
B is about 4 per cent more massive than CM Dra 1, but their 
measured radii are the same within the error. The mass- 

1 available online at ftp://phoenix.hs.uni-hamburg.de/NextGen/. 
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Figure 1. Observed masses and radii for the FC stars (show- 
ing the error bars) in the two systems compared with theoretical 
mass— radius tracks at 4 Gyr, calculated with different metallici- 
ties and a = 2.0. 



radius relationship followed by CM Dra 1 and 2 is indeed 
very different from that of the FC components of KOI 126, 
as Fig. [T]shows. Correspondingly, the effective temperatures 
of CM Dra are significantly lower than the range estimated 
for KOI-126 B and C. From the raw observational data, 
taken at face value, we can expect the best-fitting models 
for the FC stars in the two systems to have quite different 
properties. 

In modelling the two systems, we used different strate- 
gies to cope with their different characteristics and available 
observational constraints (e.g., spectroscopy of KOI-126 A, 
T cfl of CM Dra 1 and 2, etc.). For KOI-126, we looked for a 
self-consistent solution (in age and composition) for all the 
three components, without making any a priori assumption 
on the age, as is explained in detail in Section [4. II For CM 
Dra, on the other hand, we adopted the age estimate from 
the literature, since the age sensitivity of FC models is too 
weak to provide a meaningful self-consistency check (e.g., 
the radius variation is of the order of 2-3 per cent over the 
whole main sequence lifetime, 1-10 Gyr, and about 1 per 
cent between 1 and 5 Gyr). 



3 CONSTRUCTION OF THE STELLAR 
MODELS 

All the models were calculated us ing the Yale Rotating S tel- 
lar Evolution Code (YREC, see iDemarque et al.l 20081) in 
its no n-rotating configuration. We used the Ferguson et al.l 
|2005l ) opacities at low temperatu res, the OPAL Rosselan d 
opacities at high temperatures (jlgles ias & Ro gers! 19961), 
and t he OPAL 2005 equation of state ( Rogers fc Navfonovl 
120021 ). The energy generati on rates are calculated a ccord - 
ing to the prescription of iBahcall fc PinsonneauhJ (| 19921 ) 
and the diffusion co efficients for h elium and heavy elements 
are those given by iThoul et al.l Jl994h. Convection is de- 
scri bed using the MLT (jBohm- Vitensd fl95cl ) . We adopt 
the iGrevesse fc Sauv al ( 1998) value of the solar metallicity, 
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Figure 2. Upper panel: P-T atmospheric profile for a 0.2 Mq, 
4 Gyr-old star, calculated with Y = 0.28, [M/H] = 0.0, and a = 
2.0. Different T-r relations are used in the atmosphere. Lower 
panel: the same as in the upper panel, but with [M/H] = —0.5. 



The integration of atmospheric layers in YREC requires 
that the temperature -optical depth (T-r her eafter) relation 
is specified. While the lKrishna Swamvl (1 19661 ) (KS in the fol- 
lowing) or Eddington grey approximations give reasonably 
accurate results for the Sun or more massive stars, respec- 
tively, this is not true for the FC models (see section 2.5 of 
IChabrier fc Baraffei Il997l . and also Section O below) . For 
this reason, in all the modelling discussed here, we used the 
T-r relations taken from the NG model atmospheres of ap- 
propriate (logfl,T B fi). 

As a reference, with the input physics described so far, 
a standard solar model calculated with YREC has an initial 
helium content of Yn = 0.278 and a value of the MLT pa- 
rameter of a = 1.82 (a — 2.14) when using the Eddington 
(KS) T-r relation. 



3.1 The impact of the T-r relation on FC models 

The global parameters of a FC star (i.e. R, L, T e ff) are highly 
sensitive to the treatment of the outermost layers. In partic- 
ular, the model atmosphere used to determine the surface 
bound ary conditions plays a cru cial role. As was pointed 
out by IChabrier fc Baraffei (| 19971 ). a T-r relation based on 
a grey atmosphere and/or radiative equilibrium is incorrect 
for Teg < 5000 K, because the recombination of molecu- 
lar hydrogen in the envelope reduces the entropy and the 



4 F. Spada and P. Demarque 



adiabatic gradient, thus favouring the penetration of con- 
vection in optically thin layers. Moreover, a grey treatment 
produces cooler atmospheric profile s below the photosph ere, 
thus leading to overestimated T e g (jChabrier et al.|[l996T ). 

We investigated the impact of these effects on our calcu- 
lations by comparing models of FC stars constructed using 
the Eddington, KS and NG T-r relations. Fig. [5] shows the 
pressure-temperature profile in the outer layers of a star of 
0.2 Mg with so lar and moderately metal -poor composition 
(cf. fig. 5a-b of IChabrier fc Baraffd Il997r i . The models cal- 
culated with the grey T-r relations have a discontinuity in 
the temperature gradient near the photosphere (i.e. around 
T — T c ff). This is due to the abrupt truncation of convection 
at the transition from the envelope to the atmosphere, where 
radiative equilibrium is artificially enforced throughout. The 
models using the NG temperature profile, on the other hand, 
have a much shallower temperature gradient and a smooth 
envelope-atmosphere transition. The net effect is that, using 
Eddington (KS) grey atmosphere models, we obtain effec- 
tive temperatures which are overestimated of about 200 K 
(fOO K) with respect to the NG models. 

Convection is treated according to the standard MLT 
in the PHOENIX ID code. Since the models in the NG 
grid were all constructed with a = f.0, we were unable 
to estimate directly the impact of the value of a in the 
atmosphere on the global stellar quantities. Moreover, a 
value of a = 1.0 in the atmosphere is larger than the one 
needed to achieve a satisfactory fit of s pectral line profiles 
(Ivan't Veer-Menneret fc Megessierlll996l . see also section 2 of 
IPiau et al.l 2011 ). In fact, inconsistent values of a between 
the atmosphere and the interior of the same stellar model 
are often required to match satisfactorily all the available ob- 
servational constraints. Numerical simulations of convection 
also show that assuming the same value of the mixing length 
throughout all the o uter stellar layers m ight be an oversim- 
plification (see, e.g. iTanner et al. 1201 ll ). In our modelling, 
however, this inconsistency does not seem to affect signifi- 
cantly the temperature profile, as is apparent from Fig. [2] 



4 RESULTS 

4.1 Modelling the system KOI-126 

For the system KOI-126, we tried to achieve a self- 
consistency in the best-fitting models of the three compo- 
nents, as follows. Given its mass and the pr eliminary es- 
timate of the age of the system reported by ICarter et al.l 
(|201ll ) of 4± 1 Gyr, KOI-126 A is probably in the early post- 
main sequence phase, on the way to the red giant branch. 
Helium and heavy element diffusion and core overshooting 
were taken into account when modelling this star. With fixed 
values of the core overshooting parameter a ov and metal 
content [M/H], we adjusted the helium content and the 
MLT parameter in order to produce the model of KOI-126 A 
with the closest agreement with the observed T e ff,A and Ra 
(we choose to use the T e ff-R plane to take advantage of the 
high precision in the radius determination). The resulting 
best-fitting model of the primary star gives the most strin- 
gent constraint on the age of the system, since the low-mass 
components undergo very small changes once they are on 
the main sequence. We then tried to fit the radii of KOI-126 



Table 2. Best-fitting parameters for KOI-126 A (ages in Gyr). 
The meaning of the a and age ranges is explained in the text. 









[M/H] 


= +0.07 








a ov 


0.05 


Q ov 


= 0.10 


Q ov = 


0.20 


Y 


a 


age 


a 


age 


a 


age 


0.25 


1.0—1.3 


3.6—3.8 


1.1—1.3 


3.9—4.0 


1.3—1.4 


4.0—4.2 


0.28 


1.0 


2.8 


1.1-1.3 


3.0-3.1 


1.1-1.3 


3.0-3.2 


0.30 


< 1.0 


2.3 


< 1.2 


2.5 


<1.3 


2.6 








[M/H] 


= +0.15 










0.05 


Ofov 


= 0.10 




0.20 


Y 


: 


age 


a 


age 




age 


0.25 


1.5-2.0 


4.4-4.7 


1.2-2.0 


4.2-4.7 


1.5-2.0 


4.7-4.9 


0.28 


1.0-1.1 


3.1-3.2 


1.2 


3.4 


1.3-1.4 


3.6 3.7 


0.30 


1.0-1.2 


2.6-2.8 


1.0-1.2 


2.6-2.8 


1.0-1.5 


2.6-3.2 








[M/H] 


= +0.23 










0.05 


<%ov 


= 0.10 




0.20 


Y 


: 


age 


a 


age 




age 


0.25 


1.7-2.0 


5.2-5.3 


1.7-2.0 


5.1-5.3 


1.5-2.0 


5.1-5.4 


0.28 


1.2-1.3 


4.1-4.2 


1.2-1.5 


3.7-3.8 


1.3-1.5 


3.9-4.0 


0.30 


1.2 


3.2 


1.2-1.3 


3.2 


1.4 


3.4 



B, C for each composition, determining an age range that 
can be used to check for compatibility with the model of the 
primary. 

In modelling KOI-126 A, which is massive enough to 
have a convective core, we tried different values of the core 
overshooting parameter, within the range a ov = 0.05-0.2. 
This range is typically a dopted for early post-main sequence 
subgiant stars (see, e.g., | Demarque et al.ll2004l . IClaretJl2007l . 



iDeheuvels fc Michell 12011 ). Interestingly, for this star, we 
found a very close agreement between the evolutionary 
tracks calculated with NG (T off = 5800 K, logp = 4.0) and 
Eddington (grey) T-r relations. 

Table[2]reports the parameters for the best-fitting mod- 
els. With the values of q ov , [M/H], and Y held fixed, the 
MLT parameter a determines the evolutionary path of the 
star in the (T e s, R) plane. The range of a shown corresponds 
to evolutionary tracks that cross the observational box in 
the (T e s,R) plane, i.e., it corresponds to the models that 
passed within 1 a or less from both Ra and T e fr,A- Since the 
radius determination is much more precise compared with 
that of the T e g, our age estimate is that of the model which 
has R = Ra- These intervals of a and age do not have the 
meaning of formal error bars, but are reported here to give 
a quantitative idea of the sensitivity of the models to the 
various parameters. 

In general, a ov influences the best-fitting value of a and 
therefore has an indirect impact on the age (lower a values 
result in younger age estimates). As for the composition, 
higher helium content produces younger models, while a 
higher metallicity requires larger values of a and leads to 
an increase in the age of the models. 

We calculated models of the FC components of KOI-126 
with each possible choice of (Y, [M/H]) shown in Tabled 
keeping the MLT parameter fixed with the value a = 2.1. We 
did not make any attempt to account for the T e g constraint, 
since only a very rough estimate is available for these stars. 
Table [3] shows the age interval compatible with the observed 
radius for each composition. Consistent solutions (marked in 
bold in Tables [2] and [3]) exist for the following compositions: 
[M/H] = +0.15, Y = 0.28 and Y = 0.30; [M/H] = 0.23, 
Y = 0.25; the ages are 3.4 ± 0.3 Gyr, 2.8 ± 0.4 Gyr and 
5.2 ± 0.2 Gyr, respectively. In conclusion, using the central 
determination of [M/H] and a value of Y = 0.28, close to 
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Table 3. Age intervals compatible with the observational con- 
straints on the radius for KOI- 126 B and C. All the models were 
calculated with a = 2.1. 





\M/H] 


= +0.07 


Y 


KOI-126 B 


KOI-126 C 


0.25 


> 4.5 




0.28 


> 3.0 


> 4.2 


0.30 


2.2-4.5 


> 3.0 




[M/H] 


= +0.15 


Y 


KOI-126 B 


KOI-126 C 


0.25 


> 3.5 


> 5.0 


0.28 


2.0-4.5 


> 2.7 


0.30 


< 3.5 


1.5-5.0 




[M/H] 


= +0.23 


Y 


KOI-126 B 


KOI-126 C 


0.25 


2.0-5.5 


> 2.7 


0.28 


< 3.7 


1.5-5.5 


0.30 


< 3.0 


< 4.0 



the solar calibration, we obtain a consistent solution for the 
system, with an estimated age range of 3.1-3.7 Gyr. 

4.2 Modelling the system CM Dra 

For CM Dra, we assumed the age quoted in the literature, of 
about 4 Gyr, and we allowed for a variation of the chemical 
composition in a very broad range, to take into account the 
high level of uncertainty in its determination. 

As was already clear from Fig. \T\ it is impossible to rec- 
oncile the models of CM Dra with the observed values of 
the radii when using a metallicity in the range quoted in the 
literature ( — 1.0 < [M/H] < —0.6, see also the upper panel 
of Fig. [3}. In fact, even allowing all the parameters to freely 
vary within reasonable limits, we always find that our theo- 
retical models underestimate the radii (and overestimate the 
T c ff). The upper panel of Fig. [3] also shows that, to bring the 
theoretical models almost within 1 a of agreement with the 
observations, a rather extreme metal-rich composition must 
be used {[M/H] = +0.5, or Z = 3 Z©). A moderate ef- 
fect of the helium content is also apparent from the upper 
panel of Fig. [3] with higher values of Y resulting in slightly 
higher effective temperatures and larger radii. Such a high 
metallicity is unlikely for a system of Population II, but it 
should be emphasized that this membership attribution is 
based only on the circumstantial evidence provided by the 
proper motion and that the direct estimate of the metallic- 
ity is very controve rsial (see the discussion in section 5 of 
iMorales et"al]|2009h . 

The possibility of non-standard effects in the structure 
of CM Dra components has been invoked in the past to ex- 
plain the radius anomaly. Both stars are usually assumed to 
have synchronised their spin with the orbital period, which 
makes of them quite fast rotator s. They ar e also known 
to have spotted surfaces (see, e.g. lLacvl[l977l ). which is an 
indicator of the presence of subsurface magnetic fields. A 
fast rotation and/or highly magnetised re gime can reduce 
the efficiency of th e convective transport ( Goueh fc Tavlen 
19661; lTavlerill973h. This effect was taken into account by 



Chabrier et al.l (|2007m . as well as a simplified parametriza- 



tion of the star-spots. They found that both effects can be 
phenomenologically represented by a reduction of the MLT 
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Figure 3. Upper panel: comparison of the observational data 
for CM Dra with theoretical M-R relations at 4 Gyr calculated 
with different compositions and a = 2.0. Lower panel: as above, 
but keeping the composition fixed at solar values (i.e. Y = 0.28, 
Z = Zq) and using different values of a. 



parameter a. We thus tried to construct models of CM Dra 
with reduced a, which suffices to the purpose of investigating 
such activity-induced effects within the standard theoretical 
framework, keeping the number of free parameters to a min- 
imum. As the lower panel of Fig.[3]shows, with Y = 0.28 and 
solar metallicity, the disagreement can be brought within 
about 2 cr if a value of a as low as 0.2 is used. However, 
even reducing a down to 0.05, we were unable to construct 
models with stellar radii compatible with the observations 
within the 1 a uncertainty. For a < 0.05, the MLT treatment 
implemented in YREC gives rise to numerical instabilities. 
In this sense, tuning a as a free parameter is insufficient, 
by itself, to reconcile the models of CM Dra 1 and 2 with 
the observations. As the moderate sensitivity of FC stars 
to a requires such a drastic, possibly unphysical reduction, 
we also tried to model CM Dra using a combination of the 
composition (Y and Z) variation and a reduction effects. 
With [M/H] = +0.3 (i.e. Z = 2Z e ) and Y = 0.28-0.30, for 
example, we find that the discrepancy is lower than 2 a for 
a — 1.0. Clearly, a better determination of the composition 
of the system, based on more robust evidences than kine- 
matic arguments, is required to conclusively settle this issue. 
Since this information is currently missing, we refrain from 
attempting a precise modelling of CM Dra, as was done in 
Section [4.11 for KOI-126. We believe, nonetheless, that our 
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Figure 4. Evolutionary tracks for the radius of KOI-126 A com- 
pared with the observed range (thin horizontal lines) . The result- 
ing ages are 3.2, 4.0 and 5.0 Gyr. 



results suggest a revision of the metallicity of this system 
towards at least solar-like (if not higher) values. 



4.3 Comparison with previous work 

Dedicated models for the system KOI- 126 were constructed , 
independently of the present work, bv lFeiden et al.l (|201lf ). 
wh o used the Dartmouth Stellar Evolution Program (DSEP, 
see lChabover et "al"]|200ll and references therein). They were 
able to match very well the observational constraints and to 
reproduce satisfactorily the mass-radius relationship for the 
FC components, using a metallicity in the quoted range and 
a solar-calibrated value of the MLT parameter. Their esti- 
mate for the age of the system is 4.1 ± 0.6 Gyr. The codes 
DSEP and YREC have a distant common origin; they dif- 
fer mainly in the treatment of some pieces of input physics, 
namely, the equation of state, the atmospheric boundary 
conditions, and the presence of turbulent diffusion. As was 
shown in the previous section, both codes are quite success- 
ful in modelling the system KOI-126; the differences between 
the two codes manifest themselves in the estimated age. To 
ease the comparison, Fig. [3] shows the evolutionary tracks 
for the radius of KOI-126 A, constructed with YREC us- 
ing the solar-calibrated MLT parameter and helium content 
(Y = 0.278, a = 1.82, Eddington T-r relation, a ov = 0.2), 
and the sam e choic e of mass and metallicity as in fig. 1 of 
iFeiden et ail (|201ll V In this way, we arrive at the age range 
4.0 ± 1 Gyr, with a sl ightly larger uncerta inty, but compati- 
ble with the results of lFeiden et al.l (|201ll h Interestingly, our 
tracks in Fig. [4] intersect the observed radius when the star 
is already out of the main sequenc e phase, whil e this is not 
the case for the models in fig. 1 of lFeiden et al.l (|201ll 'l. 



As Fig. \T\ shows, for KOI-126 B and C the theoret- 
ical mass-radius tracks at 4 Gyr in the 0.20-0.25 mass 
range, with a metallicity between solar and [M/H] — +0.30 
(Y = 0.28 and a — 2.0), match quite well the observational 
data, consistently with the detailed modelling discussed in 
Section |4~T1 

For the system CM Dra the situation is quite different. 
We find significant discrepancies in radii and T e ff when using 
q~2 and a metal-poor, solar or even moderately metal-rich 
composition (see the upper panel of Fig. [3J • Dedicated mod- 
elling of CM Dra was atte mpted in the past. Early studie s 
used polytropic models (e.g . IPaczvnski &i Sienkiewiczlll984M ; 
IChabrier fc Baraffej (|l995h , using a complete stellar mod- 
elling approach, found a good agreement between their mod- 
els \ (with [MJH] = —0.5 and Y = 0.25) and the observations 
of lLacvl (| 19771 ). The more precise data availab le today, how- 
ever, c annot be reconciled with the models of iBaraffe et al.l 
ljl998h . the discrepancy being worse for a m etal-poor com- 
posit ion with respect to a solar composition ((Morales et al.l 
l2009h . Modelling non-standard input physics, such as stel- 
lar a ctivity and /or magne tic fields, has been attempted as 
well ()Chabrier et al.ll2007M . These effects, however, need not 
be the only (nor the dominant) explanation. The composi- 
tion, in fact, which is very poorl y constrained, can h ave a 
significant impact on R and T e ff ijBurrows et al.ll201ll ). 

In conclusion, our models match satisfactorily the ob- 
servational data for all the components of KOI-126. With 
[M/H] = +0.15 and Y = 0.28, we obtain an age of about 3.6 
Gyr (based mainly on the modelling of KOI-126 A), which is 
compatible with the i nitial estimates of Carter et al.l (|201 lh 
and with the results of lFeiden et al.l l|201ll ) within the errors. 
We are unable, on the other hand, to reconcile our models 
with the observed radii and T c s of CM Dra if we assume the 
currently estimated metal-poor composition. We find that 
the MLT parameter has a modest effect on the structure of 
FC stars, which is nevertheless significant with respect to 
the very precise measurements available for these systems 
for a < 1.0. For example, the relative change in radius for a 
star of 0.24 M© is up to 2.5 per cent when a goes from 2.0 
to 0.1. The agreement of CM Dra can therefore be improved 
using both higher metallicity and a reduced a (e.g. it can 
be brought within 2 a for Z — Zq and a — 0.2 or Z = 2 Zq 
and a = 1.0). 
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5 DISCUSSION AND CONCLUSIONS 

We modelled the systems KOI-126 and CM Dra with the 
YREC code, using up-to-date input physics. Both systems 
have very precise determination of stellar masses and radii. 
This allows to test in detail the theoretical models of the FC 
components, which is the main goal of this work. 
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